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Abstract. A census of massive galaxies at redshift increasingly higher than z ~ 1 may provide strong constraints 
on the history of mass assembly and of star formation. Here we report the analysis of three galaxies selected in 
the Hubble Deep Field South at Ks< 22 on the basis of their unusually red near-IR color J-K> 3. We have used 
population synthesis models to constrain their redshifts and their stellar masses. One galaxy (HDFS-1269) is at 
redshift z phot ~ 2.4 while the other two (HDFS-822 and HDFS-850) are at z phot ~ 2.9 - 3.0. All three galaxies 
have already assembled a stellar mass of about 10 11 Mq at the observed redshift placing the possible merging 
event of their formation at z>3.5. The inferred mass weighted age of their stellar populations implies that the 
bulk of the stars formed at zf > 3.5. The resulting co-moving density of Mstars^W 1 Mq galaxies at (z) ~ 2.7 is 
p = 1.2±0.7xl0 -4 Mpc -3 , about a factor two higher than the predictions of hierarchical models. The comparison 
with the local density of galaxies implies that the three galaxies must have already formed most of their stellar 
mass and that they cannot follow an evolution significantly different from a passive aging. The comparison with 
the density of local L>L* early types (passively evolved galaxies) suggests that their co-moving density cannot 
decrease by more than a factor 2.5-3 from z = to z ~ 3 suggesting that up to 40% of the stellar mass content 
of bright (L>L*) local early type galaxies was already in place at z > 2.5. 

Key words. Galaxies: evolution; Galaxies: elliptical and lenticular, cD; Galaxies: formation 



1. Introduction 

Deep near-IR surveys are unveiling sources with unusually 
red near-IR colors (J-K> 3). They are extremely rare at 
magnitudes brighter than K=20 while their surface den- 
sity increases at fainter magnitudes. Only one source red- 
der than J-K=3 is present in the the 65 arcmin 2 surveyed 
by Hall et al. (2001) down to K~19.5 while 5 of them 
appear at Ks< 21 over a sub- area of 43 arcmin 2 of the 
ESO Imaging Survey (Scodeggio and Silva 2000). One red 
object (HDFN-JD1) was found in the Hubble Deep Field 
North (HDFN; Dickinson et al. 2000). It has a magnitude 
K~ 22 and no counterpart at wavelengths shorter than 1.2 
/im. Maihara et al. (2001) and Totani et al. (2001) noticed 
the presence of four sources at magnitudes K'>21 with 
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color J-K> 3 in the Subaru Deep Field (SDF). Objects 
with these unusually red near-IR colors were also noticed 
in the HDF-S by Saracco et al. (2001) at K> 20.5. 

The nature of these sources has not yet been firmly 
established even if it is quite certain that they are not 
galactic objects. Indeed, very low mass stars, such as L- 
dwarfs, can display colors only slightly redder than J-K=2 
(Chabrier et al. 2000; Kirkpatrick et al. 2000). Stars heav- 
ily reddened by circumstellar dust due to undergoing mass 
loss, such as Mira variables and carbon stars, can be red- 
der than L-dwarfs. However, Whitelock et al. (1995; 2000) 
found only 2 stars having J-K~ 3 out of the 350 Mira 
and mass-losing stars observed. It seems unlikely that the 
extremely small fields of the HDFs and of the SDF can 
contain so many extremely rare stars at high galactic lat- 
itude. Furthermore, the apparent K-band magnitude of 
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these unusually red objects (K>20) would place them out 
of the Galaxy at a distance larger than 5 Mpc if they were 
stars. 

Cutri et al. (2001) and Smith ct al. (2002) find 4 QSOs 
with colors J-K> 3 among the 231 red AGNs selected 
using a J-K> 2 criterion from the Two Micron All-Sky 
Survey (2MASS). They are brighter than K=14 and are 
at z < 0.3. If the unusually red sources seen at K>20 were 
dominated by AGNs at these redshifts they would be at 
least 10 3 times less luminous then the 2MASS AGNs, i.e. 
too faint to be AGNs. On the other hand, un-obscured 
AGNs at larger z would get rapidly bluer since the rest- 
frame near-IR excess would be redshifted beyond the K- 
band. In the case of dust obscured AGNs, they should be 
reddened by at least Ay ~ 3 mag and placed at z > 2 
to match the observed J-K color. However, such values of 
extinction characterize AGNs for which the rest-frame op- 
tical luminosities are usually dominated by the continuum 
of the host galaxy (e.g. Maiolino et al. 2000). 

The Extremely Red Objects (EROs) studied so far 
(e.g. Thompson et al. 1999; Cimatti et al. 1999, 2002; 
Daddi et al. 2000; McCarthy et al. 2001; Martini et al. 
2001; Mannucci et al. 2002; Miyazaki et al. 2002) are char- 
acterized by near-IR colors usually bluer than J-K~ 2.5. 
Colors redder than J-K^ 2.5 arc not expected even for 
passively evolved galaxies down to z ~ 2 (e.g. Saracco ct 
al. 1999). Indeed, all the EROs spectroscopically observed 
so far lie at z < 2 (e.g. Cimatti et al. 2002; Saracco et al. 
2003). Thus, the unusually red near-IR color character- 
izing these objects suggests redshifts z>2 and, possibly, 
a component of dust absorption. Dickinson et al. (2000) 
consider various hypothesis for the nature of HDFN-JD1: 
from the most extreme of an objects at z>10, justified 
by the non-detection of the object from 0.3 to 1.1 /im, to 
the least extreme of a dusty galaxy at z > 2 or a max- 
imally old elliptical galaxy at z>3. The analysis of Hall 
et al. (2001) suggests a redshift z ~ 2.4 for their unusu- 
ally red object. Totani et al. (2001), by comparing the J-K 
color and the surface number density of the red sources 
in the SDF with model predictions, conclude that they 
are best explained by dusty elliptical galaxies at z ~ 3 in 
the starburst phase of their formation. Thus, the analy- 
sis of unusually red near-IR objects performed so far place 
these galaxies at z > 2 — 3. Our knowledge of the Universe 
at these redshifts comes mostly from the Lyman-Break 
Galaxies (LBGs) selected through the U-dropout method 
based on UV-optical color (e.g. Steidel et al. 1996). The 
unusually red objects above are missed by this selection 
technique both due to their faintness at optical wavelength 
and to their different optical colors (see e.g. Vanzella et al. 
2001). Consequently, also the information they bring rele- 
vant to the Universe at that z are missed. For this reason 
they could be extremely important to probe the Universe 
at these redshifts. 

In this paper we present the analysis based on a multi- 
band data set (from 0.3 /im to 2.15 /im) of three J-K> 3 
sources selected at Ks< 22 on the HDF-S. The near-IR 
data have been collected by the Faint Infra-Red Extra- 



galactic Survey (FIRES, Franx et al. 2000). In §2 we 
present the imaging, the photometry and the spatial ex- 
tent analysis of the three sources. In §3 we derive the red- 
shift and, consequently, the stellar masses of the galaxies 
through the comparison of the data with population syn- 
thesis models. In §4 we derive the co-moving spatial den- 
sity of these objects and we try to constrain their forma- 
tion and evolution in §5. We summarize our results in §6. 
Throughout this paper, magnitudes are expressed in the 
Vega system unless explicitly stated otherwise. We adopt 
an Cl m = 0.3, CIa = 0.7 cosmology with H o =70 km s" 1 
Mpc" 1 . 

2. Imaging and photometric analysis 

The optical images are the Version 2 of the Hubble 
Space Telescope images in the F300W, F450W, F606W 
and F814W filters (U300, B 450 , V@o6 and L314 hereafter, 
Casertano et al. 2000). Near-IR images centered in the 
HDF-S come from the FIRES project (Franx et al. 2000) 
and they have been obtained with the ISAAC spectro- 
imager at VLT-Antu in the filters Js, H and Ks. The three 
sources analyzed here were originally noticed by Saracco 
et al. (2001) in the analysis of a first set of such data (~ 8 
hours of exposure) and turned out to be EROs (I-K> 4) 
from the analysis of the multi-band data by Vanzella et 
al. (2001; object IDs: 822, 850 and 1269 hereafter). The 
photometry and the analysis presented here are based on 
the final FIRES data set (Labbe' et al. 2003) summing 
up to about ~ 30 hours exposure per filter with a mea- 
sured FWHM~ 0.55 arcsec. The near-IR data have been 
reduced following the recipe described in Saracco et al. 
(2001). Photometric redshifts and analysis of the whole 
multi-band data sample is presented in Fontana et al. 
(2003a). In Fig. 1, the optical (B 450 , V606> I814) and the 
near-IR (Js, H and Ks) images of the three sources are 
shown. The intensity levels in each image are constrained 
within -lcr and 3cr from the background. Colors have been 
measured within the Ks-band detection isophote by us- 
ing Sextractor (Bertin & Arnouts 1996) in double image 
mode. To this end, we first rebinned the IR images to the 
same pixel size and orientation of the WFPC2 images. 
Then, the WFPC2 images were smoothed to the same ef- 
fective PSF of the IR images as described in Vanzella ct al. 
(2001). Photometric errors have been obtained by rescal- 
ing the formal SExtractor errors by a correction factor 
to take into account the correlation of neighboring pix- 
els. The procedure used to estimate the correction factor 
is described in Vanzella et al. (2001). None of the three 
sources has a detectable flux in the Usoo-band, while only 
one of them (HDFS-1269) has a reliable flux in the B450 
band. The photometry we obtained is in good agreement 
with the one derived by Labbe' et al. (2003). In Tab. 1 
the optical and near-IR colors of the three sources are re- 
ported. It is worth noting that red AGNs have optical-IR 
colors 1.5-2 mag bluer than those of our sources (Francis 
et al. 2000; Smith et al. 2002). Thus, even if the presence 
of an AGN cannot be excluded (see e.g. van Dokkum et al. 
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Table 1. Photometry of the three EROs 1269, 822 and 850. Magnitudes arc in the Vega system. Colors are measured 
within the aperture defined by the Ks-band isophotc. The Ks-band magnitude is the MAG_AUTO measured by 
Sextractor. 



Object 


RA(J2000) 


Dec(J2000) 


U 3 oo-Ks 


B450-KS 


Veoe-Ks 


I 8 i4-Ks 


Js-Ks 


H-Ks 


Ks 


HDFS.1269 


22:32:49.25 


-60:32:11.62 


>7.7 


7.6±0.5 


6.9±0.2 


5.5±0.1 


3.0±0.2 


1.2±0.1 


20.80±0.05 


HDFS.822 


22:32:53.39 


-60:31:54.33 


>7.6 


>7.6 


7.0±0.3 


5.7±0.2 


3.5±0.3 


1.4±0.2 


20.9±0.1 


HDFS_850 


22:32:53.01 


-60:33:57.03 


>6.8 


>6.8 


7.1±0.3 


5.6±0.2 


3.7±0.4 


1.1±0.2 


21.7±0.1 




Fig. 1. From left to right: B 450 , V606> 1814, Js, H and Ks band images of the three EROs 1269 (upper panel), 822 
(middle panel) and 850 (lower panel). The images are 5x5 arcsec centered on the sources. The intensity of each image 
is in the range [-lcr;+3cr] from the background. The lines shown for reference in the images of 1269 are 1 arcsec width. 



2003), the optical-to-NIR colors exclude, once more, that 
the continuum emission is dominated by nuclear activity. 



The three galaxies appear very compact in the K-band 
image. In Fig. 2 we compare the K-band surface brightness 
profile of the three galaxies with the profile of a bright 
star in the field. The two brightest galaxies (HDFS-1269 
and HDFS-822) are clearly resolved while HDFS-850 is 
not resolved. Assuming that it is at a redshift comparable 
to the other two sources, as suggested by its colors, it 
is ~1 mag intrinsically fainter. Using the magnitude-size 
relations in Snigula et al. (2002) we found that, for the 
same profile, HDFS-850 would have a scale radius ~1.6 
and -3 times narrower than HDFS-822 (or HDFS-1269) 
in the case of an exponential and a de Vaucouleur profile 
respectively. Thus, considering the angular extent of the 
two brightest galaxies, it is expected that the profile of 
HDFS-850 is dominated by the seeing. 



3. Photometric redshifts and stellar masses of 
J-K>3 galaxies 

In this section we compare the data with population syn- 
thesis models to constrain the redshift and the stellar mass 
of the three galaxies. The grid of templates used is based 
on the latest version of the Bruzual & Chariot (1993) mod- 
els and it is similar to the one described in Fontana et 
al. (2003a, b) in terms of star formation histories (SFHs), 
metallicity and extinction (Ay). The SFHs considered, 
besides the simple stellar population (SSP) and the con- 
stant star formation rate {est), are described by an ex- 
ponentially declining star formation rate (SFR) with c- 
folding time (r) in the range 0.1-15 Gyr. In this paper, 
besides the set of synthetic templates obtained with the 
Salpeter (Sal; Salpeter 1955) initial mass function (IMF; 
0.1M© < M < 100M Q ), we considered also three analo- 
gous sets generated with the Miller-Scalo (MS; Miller & 
Scalo 1979), Scab (Sea; Scalo 1986) and Kroupa (Kro; 
Kroupa 2001) IMF respectively and a set based on the 
Kro IMF derived from the spectrophotomctric models of 
Chariot & Longhetti (2001). These latter take into ac- 
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Fig. 2. Radial surface brightness profile of the three EROs 
discussed here (points) compared to one isolated star in 
the HDF-S (line). The profiles are normalized to the flux 
within the first aperture. 

count consistently the emission from stars and from the 
gas (Kro+em) and allow to probe the possibility that the 
extreme near-IR colors are due to emission lines rather 
than to continuum emission (see e.g. Totani et al. 2001 
and Franx et al. 2003). The parameters used to generate 
the grid of synthetic templates are summarized in Tab. 2. 

The x 2 -minimization procedure of hyperz (Bolzonella 
et al. 2000) has been applied to find the best fitting spec- 
tral template to the observed colors for each set of IMF, 
SFH and metallicity. In practice, for a given IMF, metal- 
licity and SFH, we find the best fitting template, i.e. the 
corresponding redshift and age. In the best fitting proce- 
dure the extinction has been allowed to vary within the 
range < Ay < 2 mag and, at each z, galaxies have been 
forced to have age lower than the Hubble time at that z. 
We repeated this procedure for the 5 IMFs, the 3 metal- 
licity values and the 10 SFHs considered. Among the 150 
best fitting templates thus obtained we considered those 
having a probability P(x 2 )> 0.68 (where x 2 is the reduced 
X 2 ). For each galaxy, this set of templates has been used to 
define the range of variability of the relevant parameters 
we derived: z, SFHs and relevant age, stellar mass. 

For each template considered, we derived the stellar 
mass M stars of the galaxy from the mass-to-Ks-band light 
ratio since it is relatively insensitive to the star formation 
history with respect to the optical bands (Chariot 1996). 
The mass in this ratio (and thus M s tars) results, in fact, 
by the integral of the SFR over the SFH of the template 
in the interval < t < t age , where t age is the age of 
the template. This mass corresponds to the processed gas 



Table 2. Parameters used to construct the grid of tem- 
plates 



SFH r [Gyr] 


0.1, 0.3, 0.6, 1, 2, 3, 5, 15, SSP, est 


Metallicity 


0.2 Z 0> 0.4 Z 0) Z Q 


A_V [mag] 


0^2 


Extinction law 


Calzetti et al. (2000) 


IMFs 


Sal, Sea, MS, Kro, Kro+em 



involved in the star forming process in this interval (see 
e.g. Madau et al. 1998). Thus, this is the mass that has 
been burned into stars at some time in this interval, i.e. the 
mass that was contained and/or is still locked into stars. 
We have also considered the mass locked into stars (M.*) 
at the age of the galaxy as resulting from the difference of 
■Mstars and the processed gas returned to the interstellar 
medium through stellar winds. As far as our three galaxies, 
we have found that M* > 0.8M stars for Sal IMF and 
M* > 0.6M s tars for MS IMF, this latter being the lowest 
fraction of mass locked into stars for the same value of 
■Mstars- The mass A4», is strongly related to the different 
mixture of stars and to their evolutionary stage. Thus, 
it is much more dependent on the IMF, the metallicity, 
the SFH and the extinction of the best fitting template 
with respect to M stars- For these reasons, we adopted 
■Mstars as the stellar mass of the galaxies since it is a 
robust estimate. 

The results of the best fitting procedure relevant to 
the Sal, MS and Kro IMFs are summarized in Tab. 3. We 
do not report the results obtained with the Sea IMF since 
the best fits obtained with this IMF were always worst 
than those obtained with the other IMFs and the values 
of the fitted parameters were within the ranges defined 
by the other models. The results obtained with the set 
of templates including emission lines (Kro+em) are not 
listed in the table since the best fitting templates were 
always the same as those without emission lines (Kro). 
This latter result shows that the observed extremely red 
near-IR colors are not likely to be dominated by emission 
lines. In Tab. 4, we report for each galaxy the range of 
variability of the stellar mass (AM-l^^) and of the mass- 
to-Ks-band light ratio (AM/L™ f ) defined by the fitting 
models with different IMF. 

In the following discussion of the individual objects, 
we compare the properties of local galaxies with those of 
our three galaxies. To this end, we derived a lower limit to 
their Ks-band luminosities at z = applying an evolution- 
ary correction factor. The adopted correction factor for 
each galaxy is an upper limit to the expected luminosity 
evolution in the relevant range of redshift and corresponds 
to the passive aging of the youngest stellar population 
that could populate the galaxy at the relevant redshift. 
The age considered for the stellar population is the mass 
weighted age described in §5.1. The youngest age of each 
galaxy is chosen among the relevant best fitting templates 
(P(x 2 > 0.68)). The passive aging is then traced from 
z P hot to z — assuming a SSP with age at z p hot equal to 
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0.5 1 1.5 
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Fig. 3. The best fitting templates (see Tab. 3) are super- 
imposed on the observed photometric data (filled points) 
in the B450, V606, Isuj Js> H and Ks bands of the ERO 
1269 (upper panel), 822 (middle panel) and 850 (lower 
panel) respectively. Magnitudes are in the AB system. 
Filled triangles represent upper limits. Open circles are 
the magnitudes derived from the best fitting template. 



the youngest (mass weighted) age. It is worth noting that, 
for the relevant range of redshift, the expected luminosity 
evolution is not significantly dependent on the assumed 
IMF (differences less than 0.2 mag). We considered the 
Ks-band to minimize the uncertainties in this extrapola- 
tion. Indeed, the evolutionary corrections in the Ks-band 
are much smaller and less dependent on the SFH than 
those in the optical bands. For instance, the evolutionary 
correction for a passive aging from z — 3 to z = in the 
case of a SSP is ~ 1.7 mag in Ks while it is ~ 2.8 mag in 
V. 

HDFS_1269 - This is the sole object for which Sextractor 
has reliably detected flux in the B45o-band. The best- 
fitting value to the redshift of HDFS-1269 is z p hot = 



2.4 



-0.3 

-0.2- 



The quoted errors are the formal uncertainties 
of the fitting procedure (68% confidence level). The rest- 
frame (k-corrected) absolute magnitude of this galaxy, as 
resulting from its distant modulus (46.44) is Mk — —25.1, 
i.e. L>3L*, where we considered MJ^ ~ —24.2 for local 
galaxies (Cole et al. 2001; Huang et al. 2002). The best 
fit is formally given by a 2 Gyr old, r = 0.3 Gyr model 
with an extinction Ay =0.3 and Z=0.4Z Q obtained with 
a Sal IMF (see Tab. 3). In Fig. 3 (upper panel) the best- 
fitting template is shown together with the photometric 
data points. For this galaxy 61 exponentially r > 0.3 Gyr 
decaying models fit the data (0 < Ay < 2) , providing Age 
in the range 1-4 Gyr and redshifts 1.9 < z < 2.7. The fit- 



ting models provide stellar masses of the order of 10 11 M Q 
independently on the IMF assumed, as shown in Tab. 4. 
The stellar mass derived by the best fitting templates with 
the Sal IMF is M sta rs = 1.2 x lO n M . 

The degeneracy with respect to the SFH, does not help 
in tracing the evolution of this galaxy from z p hot = 2.4 
to z = 0. Indeed, the stellar mass already assembled at 
z = 2.4, may represent from 35% (if r = 15 Gyr) to al- 
most 100% (if t = 0.3 Gyr) of the stellar mass the galaxy 
could form down to z = 0, depending on the SFH con- 
sidered among the 61 models. The lower limit to the Ks- 
band luminosity that this galaxy would have at z = is 
L z= o _L*, having estimated a passive evolutionary dim- 
ming of AMk — 1.3 mag from z p hot — 2.4 to z = 0. 

HDFS_822 - The best-fitting value to the redshift is 
Zphot = 3.0^Q 2 and the resulting absolute magnitude is 
M K ~ -26.0 (L>5L*). The best fit is formally given 
by a 1.7 Gyr old, r — 0.3 Gyr model with an extinc- 
tion Ay =0.9 and Z=0.2Z Q obtained with a MS IMF 
and it is shown in Fig. 3 (middle panel). For this galaxy, 
93 templates fit the data (0.4<Ay <2) providing Age in 
the range 0.2-2.5 Gyr and redshifts 2.4< z < 3.2 The 
stellar mass derived in the case of Sal IMF is Ai stars = 
3.3 x 10 11 M Q . For this galaxy the models provide a stel- 
lar mass always larger than 10 1:l Mq independently on the 
IMF assumed. Also in this case, the degeneracy with re- 
spect to the SFH implies that the stellar mass already 
formed at z ~ 3 may represent from 20% to ~ 100% of 
the stellar mass the galaxy could have at z=0. The lower 
limit to the Ks-band luminosity of this galaxy at z = is 
L z=0 >1.2L* (AM K ~ 1.7 mag). 

HDFS_850 - This is the faintest and the reddest of the 
three galaxies selected. In this case the best-fitting tem- 
plate is given by a 0.5 Gyr old, r = 0.1 Gyr model with 
an extinction Ay =1.9 and Z=Z Q obtained with a Kro 
IMF. This model provides a redshift z phot = 2. 9to%- For 
this ERO, 58 templates (all the r < 2 Gyr models and 
SSP model) fit well the data with extinction l<Ay <2, 
Age in the range 0.25-1.7 Gyr and redshifts 2.6< z < 3.1. 
The stellar mass derived is M sta rs = 1.4 x 10 n M Q (Sal 
IMF). The stellar mass already formed and assembled at 
z ~ 2.9 may represent from ~ 70% to ~ 100% of the stel- 
lar mass the galaxy could form down to z=0. The lower 
limit we derived to the Ks-band luminosity of this galaxy 
is L z= q c±L* (AMk — 1.5 mag). 

Extremely young ages with active star formation and 
no dust attenuation are included among the best fitting 
templates. Hence, the minimum permissible mass-to-light 
ratio for each galaxy has been derived with the proce- 
dure adopted. Consequently, the lower bound of the range 
of variability of the stellar mass given in Tab. 4 should 
be considered a lower limit to the mass in stars of each 
galaxy. Moreover, active star formation can partially mask 
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Table 3. Parameters resulting from the best fits of stellar population models to the photometric data for Sal, MS and 
Kro IMFs. The x 2 represent the reduced x 2 - The model which formally provides the best fit (the lowest \ 2 ) is marked 
with "*". The V-band absolute magnitude has been derived by the Ks-band apparent magnitude. The absolute 
magnitudes in the Ks and in the V bands are corrected for the extinction. 



ID 


SFH 


A 


P(y 2 ) 


Zphot 


Aee 

iv yi J 


Av 

i > i r, rr 


z 


My 

i > i o rr 

[magj 


M K 
[mag] 


a M/L K 

[1V1/1JJ0 


StClT S 

[1U 1V10J 




Sal IMF 






















HDFS 1 9RQ* 


r 


= 0.3 


0.33 


0.92 


2 4+0.3 

z -^-0.2 


2.0 


0.3 


0.4 


-22.9 


-25.1 


0.50 


1.2 


HDFS.822 


r 


= 0.3 


0.06 


0.99 


3.01S1 


1.7 


1.0 


0.2 


-24.1 


-26.1 


0.45 


3.3 


HDFS.850 


r 


= 0.1 


0.16 


0.98 


9 q+0- 2 


0.5 


1.9 


1.0 


-24.1 


-25.8 


0.30 


1.4 




MS IMF 






















HDFS_1269 


T 


= 0.3 


0.34 


0.92 


2 r+0.3 


2.0 


0.2 


0.4 


-22.9 


-25.1 


0.26 


0.6 


HDFS.822* 


T 


= 0.3 


0.04 


0.99 


o q+0.3 
O.U_ .2 


1.7 


0.9 


0.2 


-24.1 


-26.1 


0.23 


1.4 


HDFS_850 


T 


= 0.1 


0.17 


0.98 


2 q+o.2 


0.5 


1.8 


1.0 


-24.0 


-25.8 


0.14 


0.6 




Kro IMF 






















HDFS.1269 


T 


= 0.3 


0.38 


0.89 


2 4+0.3 

z -*-0.2 


2.0 


0.4 


0.4 


-22.9 


-25.1 


0.38 


0.9 


HDFS.822 


T 


= 0.3 


0.04 


0.99 


o n+0.3 

a.u_ 2 


1.7 


1.0 


0.2 


-24.1 


-26.0 


0.35 


2.1 


HDFS.850* 


T 


= 0.1 


0.16 


0.98 


9 q+o- 2 


0.5 


1.9 


1.0 


-24.0 


-25.8 


0.22 


1.0 



a We used M , K =3.4 (Allen 1973) 



Table 4. Ranges of stellar masses (AMl^f s ) and of mass-to-Ks-band light ratios (AM/L™ f ) defined by the fitting 
models with different IMF (Sal, MS and Kro). 



ID 


A A/f 6al 

"-^J 1 stars 

[10 n M Q ] 


AM/L|°' 
[M/L] 


stars 

[lO n M ] 


AM/L| 1S 
[M/L] 


A KA Kro 

L -*- ,v 1 stars 

[10 n M Q ] 


AM/Lf° 
[M/L] 


HDFS.1269 


1.2-2.1 


0.35-0.75 


0.6-1.2 


0.20-0.44 


0.9-1.5 


0.29-0.49 


HDFS.822 


2.3-5.6 


0.25-0.52 


1.1-2.6 


0.14-0.27 


1.7-4.2 


0.2-0.47 


HDFS.850 


1.0-1.4 


0.25-0.45 


0.5-1.2 


0.12-0.24 


0.8-1.1 


0.18-0.28 



the presence of a more massive and much older stellar pop- 
ulation (see e.g. Papovich et al. 2001) which could not be 
correctly accounted for using a single SFH as adopted in 
the present analysis. The results we obtained show that i) 
all three galaxies are at redshift 2 < z < 3; ii) they have 
a high stellar mass content (Mstars — 10 11 M Q ) already 
formed and assembled at the observed redshift; iii) they 
would populate the bright end (L z= o>L„) of the local lu- 
minosity function of galaxies even assuming they evolve 
passively. 

4. The density of J-K> 3 galaxies 

The co-moving spatial density of J-K> 3 galaxies and its 
statistical uncertainty have been estimated as 

p=Er- ^)-[E(t^) 2 ] 1/2 « 

^ v max i max 

Thus, for each galaxy we computed the volume 

rz max dy 

v — = tJ zi l^ dz (2) 

where uu is the solid angle subtended by the HDF-S area 
and z m ax is the maximum redshift at which each galaxy 
would be still included in the sample limited to Ks=22. 
Since, for a given template the k-correction in Ks is nearly 
constant for z > 2 — 2.5, in the derivation of the z max of 



each galaxy we considered the average k-correction de- 
rived by the best fitting templates at redshift z — z p hot- 
The lower bound in the integration (^i) should be defined 
by the J-K> 3 color cut. In Fig. 4 the evolution of the J-K 
color with redshift is shown for different SFHs and two dif- 
ferent values of the extinction. It has been assumed that 
the star formation begins at Zf = 6 with a Sal IMF. The 
figure clearly shows that the J-K> 3 color cut does not 
define a unique value of z^ given the degeneracy with re- 
spect to the SFH and to the extinction. Thus, we assumed 
as Z\ the lower limit to the redshift of the lowest redshift 
galaxy (HDFS-1269), i.e. z x = 1.9. We find that at the av- 
erage redshift (z) ~ 2.7 the co-moving density of galaxies 
brighter than Ks=22 with J-Ks> 3 is p = 1.2 ± 0.7 x 10~ 4 
Mpc~ 3 . In the derivation of z m ax we did not taken into ac- 
count the surface brightness dimming which would tend to 
lower the maximum distance reachable. Moreover, galax- 
ies with masses of about 1O 11 M at z > 2 can escape 
the J-Ks>3 threshold due to the uncertainties in the flux 
measure and to intrinsically bluer colors. Indeed, some of 
the galaxies selected by Labbe" et al. (2003) with color 
J-K> 2.3 at Ks< 22 turned out to have redshift z>2 and 
masses of the order of 10 11 M© (Fontana et al. 2003b). 
Thus, it is likely that this density is an underestimate of 
the density of Mstars^W^^Q galaxies at these redshifts. 
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2 3 4 

redshif t 

Fig. 4. Evolution of the J-K color as a function of redshift. 
The different curves refer to the different SFHs quoted 
in the figure. All the models have been obtained with 
Sal IMF assuming a star formation beginning at Zf — 6 
(Age=0) and two different values of extinction (Ay=0 up- 
per panel and Ay = 1.5 lower panel). The Age of the stellar 
population is shown on the top axis. The horizontal line 
defines the J-K=3 color cut used to select the three galax- 
ies. 

5. Constraining the evolution of massive galaxies 

In this section we attempt to constrain the possible evo- 
lutionary scenarios of the 3 massive galaxies. In §5.1 we 
concentrate on the formation and evolution of the stellar 
populations, while in §5.2 we consider the evolution of the 
number density of massive galaxies. 

5.1. Constraining the star formation history 

Spectrophotometric models, as those adopted in the 
present work, are not meant to describe the detailed evo- 
lutionary history of galaxies. Their aim is to describe the 
global observed photometric properties of a galaxy by 
means of an appropriate mix of stellar populations and dif- 
ferent SFHs. Real galaxies are not necessarily expected to 
follow one of these SFHs. Different combinations of them 
can provide good fits to the photometric data of a given 
galaxy implying different behavior of its past and future 
evolution (as clearly shown in §3). Keeping these limits 
in mind, we have tried to constrain the SFH of the three 
massive galaxies. 

In order to characterize the properties of their stellar 
populations we computed a more robust estimate of the 
age of the bulk of stars in the three EROs. Indeed, we 
cannot consider the age parameter provided by the best- 
fitting models given the degeneracy of the SFHs. We thus 



defined for each best fitting template a mass weighted age 
of the stars. Any template, defined by a fixed age and SFH, 
can be seen as the sum of SSPs with different ages. Each 
SSP provides the fraction of the total mass which depends 
on its own age and on the SFH describing the template 
itself. We derived the mass weighted age by summing the 
ages of the SSPs, each of them weighted on its mass frac- 
tion for each best-fitting template. The age of the bulk 
of stars in each galaxy is defined by the range spanned by 
the mass weighted age. Following this approach, we find 
that the bulk of the stars in HDFS-1269, HDFS-822 and 
HDFS-850 have ages 1.7±0.5 Gyr, 1±0.6 Gyr and 0.6±0.3 
Gyr respectively. Such estimates suggests that most of the 
stars have formed at z/>3.5 in all three galaxies. The obvi- 
ous question is how the three galaxies have formed and/or 
assembled their stellar masses, i.e. how they evolved back 
in time. 

Since in the hierarchical paradigm the age of the stars 
is not necessarily correlated with the epoch of assembly 
of the objects (Moustakas and Somerville 2002), even our 
estimate of the age of the bulk of stars does not place 
any constraint on the past assembly history (if any) of the 
three galaxies. Thus, from the properties of the stellar pop- 
ulations of the three galaxies we cannot probe their pos- 
sible formation through merging process. Anyway, their 
redshift (z ~ 3) place the possible merging event at z>3.5 
considering a dynamical time scale of 3 x 10 8 yr (Mihos & 
Hernquist 1996). On the other hand, we can probe the case 
in which the galaxies are the result of the SFHs resulting 
from the best-fitting procedure. This procedure produced 
a set of acceptable SFHs only on the basis of the photome- 
try of the three galaxies. By means of the age of the stellar 
populations, the luminosity, the stellar mass and the co- 
moving density of the three galaxies, we now try to bound 
further this set of SFHs. The following discussion and the 
results are, in fact, independent on the choice of the IMF. 
In order to obtain the ages of the bulk of stars estimated 
above in the three galaxies, the r = 0.1 Gyr model (the 
lower bound of the SFHs) requires a formation redshift 
Zf ~ 3.5. To produce a stellar mass of about 10 11 M Q 
within 0.1 Gyr (3.3<z/<3.5), an average SFR-1000 M© 
yr -1 is needed in this interval. Such high SFR resembles 
that estimated by Ivison et al. (2002) for the 36 sources 
detected with SCUBA by Scott et al. (2002) and Fox et al. 
(2002) on 260 arcmin 2 . Given the comoving density of the 
three galaxies (p = 1.2±0.7x 10 -4 Mpc -3 ) and the volume 
sampled by the HDF-S within 3.3<z/<3.5 (~ 3.1 x 10 3 
Mpc 3 ) we expect 0.4 ± 0.2 of these progenitors over an 
area equal to that of the HDF-S. This number would be 
consistent with the number of SCUBA sources expected 
on this area (<~ 0.7) on the basis of the surface density 
measured by Scott et al. Star formation time scales longer 
than 1 Gyr can in principle account as well for the stel- 
lar masses already formed and for the age of the bulk of 
stars, provided that zj > 5 and SFR< 100 M© yr -1 . On 
the other hand, Fig. 4 suggests that a time scale longer 
would not be consistent with the observed J-K color, un- 
less to invoke an extinction in excess to Ay ~ 1.5 — 2 mag. 
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However, even in this case, the star formation rate has to 
rapidly fade at z ~ 3 to match the constraints imposed by 
the local density of bright galaxies. Indeed, if the stellar 
mass already assembled (~ 10 11 M©) by the three galax- 
ies at z phot was lower than 50% of the total mass in stars 
they could form from Zf to z = 0, they would be brighter 
than 2-3L* at z = 0. Thus, at z — we should measure 
a density of galaxies brighter than 2-3L* not lower than 
1.2 x 10~ 4 . The local density of galaxies brighter than 2L* 
and 3L* is 1.8xl0" 4 Mpc" 3 and 0.5xlCT 4 Mpc" 3 respec- 
tively, i.e. comparable or lower than the density of the 
three galaxies. Thus, the photometry, the stellar masses 
already formed and assembled, the mean ages of the stel- 
lar populations together with the constraints imposed by 
the local density of bright galaxies, suggest that the three 
galaxies are consistent with galaxies which at z ~ 3 have 
formed most of their stellar mass and which will follow 
a passive aging in time soon after few Myrs, i.e. with a 
massive early type galaxy. 

5.2. Constraining the density evolution 

We have compared the density of local massive galaxies 
with the one at redshift z ~ 2.7. We have previously 
found that the luminosity of the three galaxies at z = 
has to be ^L* in the case of passive evolution and should 
not be brighter than ~2L* in the case of a slightly longer 
star formation activity. In this latter hypothesis, the three 
galaxies would represent ^60% of the whole population of 
massive local galaxies (L>2L*). In the hypothesis of pas- 
sive aging, we have to consider the local number density of 
early type galaxies. By integrating the LF of E/SO galaxies 
brighter than L* (M s tar — 10 11 M©) we obtained a den- 
sity of 3xl0~ 4 Mpc -3 , where we used <j>* E , S0 = 1.5 x 10~ 3 
Mpc -3 (Marzke ct al. 1998). Thus, the three galaxies ac- 
count for about 40% of the local population of early type 
galaxies. This result implies also that ^40% of the stellar 
mass content of local early type galaxies was already in 
place at z>2.5 and that their number density cannot de- 
crease by more than a factor 2.5-3 from z — to z ~ 3. 
In the semi-analytical rendition of hierarchical models by 
Kauffmann and Chariot (1998a,b), the predicted number 
density of M. s tar > 10 11 M© ellipticals decreases by a fac- 
tor <~ 8. Indeed, the predicted density of Ai s tar > 10 11 
M© ellipticals at z ~ 2.5 is about 5 — 6 x 10~ 5 Mpc -3 , 
a factor 2 lower than our estimate. The same conclusion 
is reached considering the predictions of Moustakas and 
Somerville (2002). They predict a number density of ha- 
los hosting local bright (L>L*) elliptical galaxies (gEs) 
at z ~ 2.5 more than a factor two lower than our esti- 
mate. The semi-analytic models of galaxy formation in 
the hierarchical clustering scenario by Baugh et al. (1998, 
2002) predict, in fact, no such massive galaxies at redshift 
z > 2. These results are summarized in Fig. 5. The figure 
suggests that the evolution of the number density of mas- 
sive (M s tars — 10 11 M©) galaxies with redshift is slower 
than that predicted by the current hierarchical models, at 



least in the redshift range < z < 3. It is worth not- 
ing that, by assuming a value of M s t ar — 10 11 M© for 
the three EROs, we derive at z ~ 2.7 a stellar mass den- 
sity p s tar = 1.2(±0.7) x 10 7 M© Mpc -3 . This density is 
<~ 30% of the total stellar mass at this redshift (Fontana 
et al. 2003b). 

The density we estimated is lower than the density 
of ellipticals (2.16±0.6 x 10~ 4 Mpc -3 ) spectroscopically 
identified at 0.85 < z < 1.3 by Cimatti et al. (2002). On 
the other hand, this estimate includes mainly ellipticals 
less massive than 10 11 M©. Indeed, given the limiting 
magnitude K=19.2 of their spectroscopic sample, they are 
sampling the counterpart of local L>0.3L* (M. s tar > 10 10 
M©) ellipticals in that redshift range. Our estimate is also 
lower than the density of Ai s tar > 10 11 M© galaxies at 
z ~ 1 estimated by Drory et al. (2001). This is expected 
since they use an approach which maximizes the stellar 
mass of the galaxies for any K-band luminosity at any 
redshift. Thus their estimate represent the upper limit of 
Mstar > 10 11 M© galaxies at that z. 

6. Discussion and Conclusions 

We have presented the analysis of the three galaxies se- 
lected on the HDF-S on the basis of their unusually red (J- 
K> 3) color. We have constrained their redshifts and esti- 
mated their stellar mass content by comparing the photo- 
metric data with population synthesis models obtained for 
different IMFs, SFHs, metallicity and extinction values. 
The grid of templates we used includes a set of template 
derived from the spectral models of Chariot & Longhetti 
(2001) which take into account consistently the emission 
from stars and from the gas. This has allowed us to ex- 
clude the possibility that the observed extreme colors are 
dominated by emission lines. We find that one of the galax- 
ies (HDFS-1269) is at z phot ~ 2.4 (P(x 2 )=0.92) while the 
other two (HDFS-822 and HDFS-850) are at z phot ~ 3 
(P(x 2 )> 0.98). All three galaxies have already assembled 
a stellar mass of about 10 11 M© at the relevant redshift. 
The fact that the two galaxies at z ~ 3 have assembled 
this mass places the possible merging event of their for- 
mation at z>3.5 considering a dynamical time scale of 
3xl0 8 yr (e.g. Mihos & Hernquist 1996). The inferred 
mass weighted age of the stellar populations places the 
formation of their bulk at Zf > 3.5 in all three galax- 
ies suggesting a substantial amount of star formation at 
these redshifts in addition to the one derived by LBGs 
(e.g. Ferguson et al. 2002). Galaxies with stellar masses 
■Mstar^W 11 M© fully assembled at z > 2 were previously 
found by other authors. For instance, Francis et al. (2001) 
find two luminous extremely red galaxies at z ~ 2.4 whose 
radial profiles suggest they are elliptical galaxies. They es- 
timate a mass ~ 10 11 M© and an age of ~ 7 x 10 8 yr for 
their stellar populations. Shaplcy ct al. (2001) find some 
galaxies with such high stellar masses and old ages among 
their LBGs at z ~ 3. Genzel et al. (2003) estimate a con- 
servative lower limit of 1.4 x 10 11 M© to the stellar mass 
of the sub-mm source at z ~ 2.8 they studied. Contrary 
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N EU M star >1011 M ( K&C 199B ) 

N gEs L>L" (M&S 2002) 




0.5 1 1.3 2 2.5 3 

z 

Fig. 5. Number density of galaxies as a function of red- 
shift. The filled circle is our estimate of the density of 
M. s tars — 10 11 M galaxies at z ~ 2.7. The filled triangle 
is the local number density of L> 2L* galaxies obtained by 
integrating the LF of Cole et al. (2001; see §4). The filled 
pentagon is the density of local L> L* early type galax- 
ies derived by the LF of Marzke et al. (1998). The open 
square is the density of L>0.3L* ellipticals from Cimatti 
et al. (2002) and the open circles are the upper limit to 
the density of M. s tars > 10 11 M Q galaxies from Drory et 
al. (2001) (see text for details). The two open triangles 
connected by the dotted line represent the range of values 
of the density of M. s tars > 10 11 M Q galaxies as derived 
by Genzel et al. (2003) from the SCUBA lens survey. The 
thin line is the density evolution of ellipticals with stellar 
mass greater than M stars — 10 11 Mq from Kauffmann 
and Chariot (1998a,b; K&C 1998); the thick line is the 
density evolution of dark matter halos hosting present- 
day bright (L>L*) ellipticals (gEs) from Moustakas and 
Somerville (2002; M&S 2002). 

to these and our findings, Dickinson et al. (2002) do not 
find galaxies more massive than 10 11 M Q at z > 2 in 
the HDFN. The small area covered by the HDFs could 
be the reason of this discrepancy even if surface densities 
of J-K>3 galaxies comparable or higher than that in the 
HDF-S are found in the other pencil beams reaching sim- 
ilar depth (e.g. Bershady et al. 1998; Totani et al. 2001). 

We found that the three galaxies must have already 
formed most of their stellar mass and that they cannot fol- 
low an evolution significantly different from a passive ag- 
ing. This suggests that J-K> 3 galaxies are most likely in 
the post starburst phase rather than in the starburst phase 
of their formation as hypothesized by Totani et al. (2001). 
These findings strongly support the thesis that J-K> 3 



galaxies are the high-z counterpart of local Mgtar^O 11 
M Q early type galaxies and agree with the recent finding 
of an increasing clustering of high-z galaxies with redder 
colors (Daddi et al. 2003; Roche et al. 2003). We estimated 
a co-moving density of galaxies brighter than Ks=22 and 
redder than J-K= 3 p = (1.2 ± 0.7) x 10~ 4 Mpc" 3 at 
the average redshift (z) ~ 2.7 which should be considered 
an underestimate of the number density of galaxies with 
■Mstar^lO 11 M at that z. This value is about a factor two 
higher than the predictions of hierarchical models rendi- 
tions by Kauffmann & Chariot (1998a,b) and Moustakas 
& Somerville (2002). In the hypothesis of passive evolu- 
tion, their luminosities at z = would be L>L*. By com- 
paring the density of local L>L* early type galaxies with 
our estimate we find that their density cannot decrease 
by more than a factor 2.5-3 from z = to z ~ 3 suggest- 
ing that up to 40% of the stellar mass contained in local 
massive galaxies was already in place at z>2.5. 
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